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Abstract 

Qh! 

(D ; It is shown that a Nambu-Goldstone (NG) field may be coherently pro- 

duced by a large number of particles in spite of the fact that the NG bosons do 
not couple to flavor conserving scalar densities like ipif). If a flavor oscillation 
process takes place the phases of the pseudo-scalar or flavor violating densi- 
' ties of different particles do not necessarily cancel each other. The NG boson 

gets a macroscopic source whenever the total (spontaneously broken) quan- 
tum number carried by the source particles suffers a net increase or decrease 
in time. If the lepton numbers are spontaneously broken such classical NG 
(majoron) fields may significantly change the neutrino oscillation processes 
in stars pushing the observational capabilities of neutrino-majoron couplings 
down to m u /300 GeV. 

PACS numbers: 14.80.Mz, 14.60.Pq 
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As is well known, spontaneous violation of global symmetries leads to the ap- 
pearance of massless Nambu-Goldstone (NG) bosons in the particle spectrum. The 
NG bosons not only have zero mass but they also have no scalar potential terms 
such as 4 . Nevertheless, they do not mediate long range forces because they only 
have derivative couplings 0. This can be understood as follows. A NG boson <ft 
associated with a quantum number A broken by the vacuum transforms under U(1)a 
as — > <p + a , a = const., and all the other fields are made invariant in the uni- 
tary gauge. The A conservation law takes the form, ignoring for simplicity possible 
mixing with other NG bosons, 

VAd^ + dpJ^O, (1) 

which is nothing but the equation of motion of the NG boson. V\ is the scale of 
symmetry breaking and is the current of the other fermion and boson particles. 
In other words, this equation reads J M = 0, where 

j» = r A + v A d^ (2) 

is an exactly conserved current. 

The source of a scalar field is of course a scalar density. The possible fermion 
bilinears are the scalar and pseudo-scalar densities like f i fj and fffsfj but they have 
to be derived from the divergence of some current J^. The actual current depends 
on the particular theory and quantum number A. is determined in leading order 
by the quantum numbers of the existing particles but receives also higher order 
corrections with a more general flavor structure. In any case, it suffices to apply the 
Dirac equation to arbitrary vector and axial-vector bilinears, written in terms of the 
fermion mass eigenstates, to conclude, as follows from 

d^Ja^lbfj = i(m + m^Jtfsfj , (3) 
<9 M Ji^fj =i(mi- rrijjfjj , (4) 

that the NG bosons couple to pseudo-scalar densities (generally these can be both 
flavor diagonal and non-diagonal) or to off- diagonal scalar densities. However cou- 
plings to flavor- diagonal scalar densities are not possible [0, [|. The pseudo-scalars 
that are diagonal in flavor vanish for free particle states and the flavor violating 
densities depend on the relative phase of distinct flavors. The natural conclusion 
has been that these relative phases cancel each other when summed over a large 
number of particles and therefore long range '1/r' NG fields are not possible. Only 
spin- dependent '1/r 3 ' interaction potentials can exist [0, [|. 
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However, there may be cases where certain quantum numbers are not conserved 
in a very large scale. For example, the total lepton number is violated if the neutrinos 
have Majorana masses, while the partial lepton numbers L e , or L T are violated 
if they are also mixed (i.e., if the non-diagonal elements of the neutrino mass ma- 
trix are non-zero in the flavor basis). In other words, this means that in physical 
processes the individual lepton numbers are violated by the neutrino oscillation phe- 
nomena, e.g. in stars or other astrophysical objects. If the unconserved quantum 
number A is spontaneously broken, it implies the existence of a massless NG boson 
- majoron [|].[] As it was pointed out by one of us [|J, in this case the non- vanishing 
divergency d^Jfc generates a classical field for the respective NG boson (majoron). 
The prototype of such system is a flux of neutrinos undergoing a flavor oscillation 
process. The key idea can be phrased as follows: as far as the neutrino oscillation 
process implies a non-conservation of the currents like J£ = v e ^^v e , etc. associated 
with the partial lepton numbers, then the conservation of the full current (P) should 
imply the existence of a classical configuration of the corresponding majoron fields 
e , etc., with non-vanishing <9 M e Q. (Clearly, a constant NG field does not carry 
much physical sense.) This seems in contradiction with what was said about the 
interaction of a NG boson with fermion scalar densities. It is the aim of this paper 
to show that it is not so because the phases of wave functions of different flavors are 
not independent from each other if transition processes take place between them. 
In that event, the wave functions of such flavors interfere constructively with each 
other and potentially form a macroscopic source of the NG boson built out of non- 
diagonal densities like fffsfj or fifj, in the basis of the mass eigenstates. To see this 
one has to be more specific about the nature of the fermion current and equations 
of motion. 

Consider a NG boson associated with the spontaneous breaking of a partial lep- 
ton number A, which could be L e , L M , L T , or any combination of these. Denoting the 
neutrinos by f a , with quantum numbers A a , and restricting ourselves to the fermion 
current, the equation of motion (|]) of the corresponding NG boson (majoron) is 
given in leading order by 

w = -i-]C^(7VAar) • ( 5 ) 

After applying the fermion equations of motion, the interactions that violate the 
lepton number A emerge in the second term. They vary from model to model 

1 We call such NG bosons majoron independently of whether the quantum number A is the total 
lepton number L or a partial lepton number L e ^. T . 
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but include in general neutrino Majorana masses. These are the lowest dimension 
terms that violate A and we confine to them here neglecting any higher order flavor 
violating gauge or Yukawa interactions. For simplicity, the neutrino mass matrix, 
m, is assumed to be real in the weak basis. Denoting by z/£, the left-handed 
neutrino fields and their charge conjugates, the equations of motion are of the form 

i^v a L = m ah vf + V^^v a L . (6) 

The potentials account for the local neutral and (Fierz-transformed) charge cur- 
rent interactions with the medium || and couplings to the NG boson namely, 

£^ = T^0^ 7 M A a z/£- (7) 



The result is 



W = Va ^ {Am)ab L ~ ul {mA)ab L 

where the quantum number A is written in matrix form. For example, one can 
consider the anomaly free quantum number L e — L^, in which case the (2-flavor) 
matrices of leptonic charge and neutrino mass have the form 

a ( 1 \ (Mi M \ . n , 

A = { -1 ) ' m= \M M 2 J ' ® 

respectively, where M is a A = L e — conserving entry and Mi , M 2 are the A- 
violating ones. The result above illustrates also that in scattering processes the 
effective coupling constants of neutrinos to majorons are essentially the neutrino 
masses divided by the scale of global symmetry breaking, V\. This point will be 
called later when discussing the observational implications. 

In order to evaluate the source terms over a system of particles one needs to 
relate the field operators with the wave functions. The chiral fields ul, are the 
left and right-handed projections of the Majorana fields 

"l + Vl = fa^j + a^ c , (10) 



where a, are annihilation and creation operators and ip single particle wave func- 
tions. The expectation value of the operators in Eq. (||]) gives then a sum over all the 
existing v particles in terms of their wave functions (75 is equal to —1 for left-handed 
spinors) : 

d^(j> = W (™A + Am) ab 75 i> h . (11) 
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The second member contains precisely the pseudo-scalar densities present in Eq. (^]). 
It remains to establish the ip equations of motion. They are non-linear at the 
operator level due to Majorana mass terms but the wave functions obey linear 



equations as follows from Eqs. dg) and (|ToD 



i^0° = m a6 V 6 -^°y*75V' a - (12) 

These are the equations of motion relevant for the most common cases of neutrino 
propagation in matter or vacuum. What is usually done is to separate the spin and 
flavor degrees of freedom by expressing the wave function as a product of a left- 
handed spinor tpQ (75^0 = — ^o)) solution of the zero mass Dirac equation 

z00o = O, (13) 

and a flavor-valued wave function (p a that obeys a well known evolution equation as 
a function of the distance travelled by each neutrino. Here, one has to go beyond 
that approximation because the pseudo-scalar densities in Eq. (|TT| ) vanish for spinors 
with well defined chirality. An approximate solution of Eq. (|1^) is 

^(*)-^ + 7 °,^ V, (14) 

which applies for the most common cases where only the scalar components Vq 
exist and even for vector potentials V a that are parallel to the neutrino velocity v. 
The wave function (p a is considered as a function of the distance travelled by the 
neutrino to the point x, s = x -v — Xq -v, while the dependence on the variable t — s 
(constant along the neutrino trajectory) is totally absorbed in ip . ip a obeys the 
familiar evolution equation |5|, |6], |7J 



ab 



where = p^/E is the neutrino velocity and V^ b = 5 a b V® . 

Using these wave functions, the majoron equation of motion (|TTD reads in leading 
order 

<p b 4tf> . (16) 



* — a t 



2 

A, — 
' 2E 



ab 

It is now clear that no flavor conserving terms contribute to the majoron source 
either in the weak basis where A is diagonal or in the mass eigenstate basis where 
m is a diagonal matrix. The existence of a majoron source relies on the interference 
between different flavor components of the neutrino wave functions. But that simply 
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means the existence of neutrino oscillations. Another necessary condition is that the 
flavor dynamics Hamiltonian does not conserve the quantum number A. Here, as a 
result of our particular specification of the neutrino interactions, A is only violated 
by the mass matrix. In general there could be also some additional non-standard 
neutrino flavor-violating interactions with matter constituents whose effects can be 
comprised in the potential in Eq. fll5|). The remarkable feature is the absence 
of cancellation due to the arbitrariness of the wave function initial phases: each 
neutrino contributes with a term that only depends on the phase invariants iplipo 
and (p a ^(p b . This crucial fact permits that in certain circumstances the contributions 
from a large number of particles add to each other with a definite sign making so a 
source of macroscopic dimensions. The precise understanding of the nature of those 
circumstances is the subject of next discussion. 

Making use of Eq. flT5| ) one can write the last equation as 

<W = ~T, *o ^V f A^ a • (17) 

Again, the second member does not depend on the initial phases of the individual 
particles. It rather depends on whether there is a net increase or decrease of the A 
number carried by the neutrinos. Take the example of a large system where only v e 
are produced out of electrons captured in nuclear reactions. The oscillations v e — > 
necessarily lead to a total decrease of the partial lepton number A = L e — that 
was initially carried by electrons. If A is a spontaneously broken quantum number, 
then the equation of motion of the corresponding NG boson 0a is (we ignore a 
possible mixing with other NG bosons |4]]): 

W A = ~e M d(i V — = - £ £ v>V ^ , (18) 

V A y ds VKv ds 

where for simplicity we have assumed that only v e <-> takes place and thus 
P Ue + P v = 1. Therefore, the shape of the majoron field a is determined in 
terms of the probability P Ve (s) = tp e > tp e /(p> (p of observing the u e flavor. Since in 
the presence of oscillations this probability is essentially smaller than P Ue {fy = 1> 
the total 0a source charge is positive. The consequence is a long range 0a classical 
field that propagates at the velocity of light. The 0a sources exist where L e — 
is either created or destroyed and the size of these regions is determined by the 
neutrino oscillation length, or resonance length in the case of resonant oscillations 
in a medium. 

An explicit example is the well known case of two flavor oscillations in vacuum |]5| . 
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With a mixing angle 9 and Am 2 mass difference the v e flavor probability is 




P^( S ) = l-sin 2 2#sin 2 —^s) , (19) 



as a function of the distance s travelled by the neutrino from its production point. 
The equation of motion of the NG field becomes 

Wa = ^ E V V sin 2 2£ ^ sin (^ a ) . (20) 

The source term changes sign as a function of s but the scale of spatial variation 
is the flavor oscillation wavelength, 4E/Am 2 , not any scale related to the number 
density or linear momentum of the particles. The neutrinos produced too far away 
give a vanishing contribution because the coherence is lost, but the ones produced 
within a sphere of the scale of the oscillation wavelength give a non-zero source term. 
Its sign depends on the distance to the point where each v e was produced but, if one 
integrates over all neutrinos, the result is finite as shows the following example. If 
the reactor or star produces v e neutrinos in a stationary basis, the NG field obeys a 
Poisson equation and the total 'charge', the volume integral of the second member, is 
positive. This can be seen from Eqs. ( jl8|) and (p!9[): the volume integration is trivial 
if one considers neutrinos propagating in plane waves and the average value of the 
probability P Ve over the neutrino spectrum converges to 1 — sin 2 26/2 at distances 
much larger than the oscillation length. For a total v luminosity equal to N u , the 
total 'charge' is 2/V\ times the number of neutrinos that undergo the transition 
v e — > per unity of time, 

xr^r^ -> "J = xr^^ 1 ~ < p -(°°)» = tf^ sin2 29 • ( 21 ) 

V A dt V A V A 

This quantity has a definite sign, regardless of the size of the reactor, the only 
assumption is the nuclear reactions feed the system with u e at a constant rate in 
time. That means that the NG field is just like the electrostatic potential due to 
an electric charge distribution: it exits at distances much larger than the size of the 
source (the place where neutrinos oscillate) with a Coulombian '1/r' shape. What 
is really different is the nature of the charge. The 'charge' of this NG field is the 
time rate of decreasing of the total L e — number carried by leptons. 

In the case of neutrino propagation in the medium, the neutrino oscillation is 
typically suppressed by the effects of the coherent interaction with the matter con- 
stituents. However, there can be a resonant neutrino conversion, so called MSW 
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(Mikheyev-Smirnov-Wolfenstein) oscillation [7|, where almost complete v e — > z/^ 
conversion takes place within the length of the resonance area. Then, the source of 
the NG field lies in the resonance shell, spheric in the case of a star, and the total 
NG charge is finite (provided that the v e — ► conversions are not compensated by 
an exactly equal number of v e — ► or — > z/ e conversions, which might happen if 
v e or ^ are also produced in the reactor or star). Consequently, a Coulombian field 
is produced with 

2 1 dN, 

Ve -> ^) (22) 



rJ1 1/ A 4vrr dt 

at radius r larger than the radius of the resonance shell. The field is constant in the 
region inside that shell. 

The fermions interact with the long range NG field through its gradient. As 
shows Eq. (0 ), d^4>\ acts as a vector potential on the neutrinos (and charged lep- 
tons). In the case illustrated above a radially moving neutrino 'feels' a potential 
energy equal to 

V NG = T^-V<p A = ±^ JL^fo - v,) , (23) 

where the upper sign is for u e , and the lower sign for u^, v e . The derivative nature 
of the NG boson couplings gives an effective vv<\> coupling constant equal to V^ 1 
over the distance to the <ft\ source. This is typically an extremely small number. 
However neutrino oscillations are sensitive to very tiny external potentials. In ad- 
dition the NG boson couplings depend on the particle quantum numbers and such 
a non-universality makes them potentially important for the neutrino oscillations 
themselves (see Eq. flI5|)). Because these are long-range fields one can conceive that 
a NG field generated in some region of a star may affect the propagation of other 
neutrinos or flavors outside that region or even out of the star. 

Some possible effects can be devised concerning the Supernova neutrinos ||. 
The NG potentials are proportional to the v luminosity over (tVa) 2 - In the case of 
Supernovae this number is comparable to the values of Am 2 /2E that are interesting 
for the solar v solutions for scales Va as large as the weak breaking scale. It was 
shown Q that the NG fields may then have a significant role in resonant v oscilla- 
tions. A special feature is the dependence of the effects on the absolute magnitude 
of the v fluxes. The signature is a surprise, i.e. oscillation patterns of Supernova 
neutrinos that are in contradiction with the solar, atmospheric and laboratory neu- 
trinos observations and/or that turn off as the v fluxes decay after the first instants 
of Supernova v emission. This may be observed with detectors capable of detecting 
Supernova neutrinos and measuring their energy spectra or time evolution. Then, 
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one could go beyond present experimental and astrophysical limits on the majoron 
effective couplings to neutrinos f| (these are typically of the order of the neutrino 
masses over the scale V\, see Eq. (§)), provided that the scale of symmetry breaking 
is under 1 TeV 0. The NG field could have some impact on solar neutrinos too 
if the scale V\ is below KeV. Though such a situation is not very plausible, there 
exist some models |TI| in which V\ can be very low and thus the effective majoron- 
neutrino coupling constants rather large, up to order 1CT 2 — 1CT 3 , still evading the 
present limits. 

The following remark is in order. The physical features of generating the classical 
majoron field due to neutrino oscillation is very much different from the mechanism 
of the majoron production due to the neutrino coherent scattering on the matter 
components of a medium - so called matter induced neutrino decay with majoron 
emission In the latter case majorons are produced as particles since the medium 



provides an energy splitting between the neutrino and antineutrino states and so the 
transitions either i/^z/ + 0orz/— >z/ + are possible, depending on the neutrino 
flavor and on the chemical content of the medium (i.e., which of two states becomes 
"heavier" in the matter background). Very low-momentum majorons could also be 
produced as a result of (stimulated) neutrino decays in the Early Universe fl2| , form- 
ing a sort of Bose condensate with large occupation numbers in the low momemtum 
part of the spectrum, but still made of individual majoron particles. On the op- 
posite, the production of the classical majoron field rather resembles the situation 
with classical electromagnetic field produced by the electric current. 

Concluding, we have shown that if a process of flavor oscillation takes place a 
pseudo-scalar density like the one of Eq. ([□]) may give rise to a long-range NG field 
as a result of the constructive interference between the wave functions of different 
mass eigenstates. Furthermore, the source of the NG boson is the time rate of 
decreasing of the quantum number associated with it, which is nothing but the term 
—dfiJ/t i n Eq. ([[]) |4j]. If the lepton numbers are spontaneously broken at the Fermi 
energy scale or below, and that means effective neutrino- majoron coupling constants 
as low as m^/300 GeV, the associated majoron fields are still significant enough for 
neutrino oscillations in Supernovae whose spectra may then show evidence of their 
existence beyond the limits || that can be reached from laboratory or astrophysical 
scattering processes. 
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